Avalanche-driven fractal flux distributions in NbN superconducting films 
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Flux distributions in thin superconducting NbN films placed in a perpendicular magnetic field have 
been studied using magneto-optical imaging. Below 5.5 K the flux penetrates in the form of abrupt 
avalanches resulting in dendritic structures. Magnetization curves in this regime exhibit extremely 
noisy behavior. Stability is restored both above a threshold temperature T* and applied field H* , 
where H* is smaller for increasing field than during descent. The dendrite size and morphology are 
strongly T dependent, and fractal analysis of the first dendrites entering into a virgin film shows 
that dendrites formed at higher T have larger fractal dimension. 



Flux jumps are known to destroy the critical state of 
type-II superconductors and suppress the apparent crit- 
ical current densityi^ In thin films the flux jumps man- 
ifest themselves in a so-called dendritic instability, i.e. 
avalanche-like penetration of magnetic flux along narrow 
branching channels. Using magneto-optical (MO) imag- 
ing the dendritic instability has been observed in super- 
conducting films of Nb)2*^ YBa2Cu307_#£i£ (triggered 
by a laser pulse) , MgB^ 8 i 9 i 10 i u and YNi 2 B 2 CJ^ Recently, 
flux dendrites were found also in films of NbaSnji 3 . a su- 
perconductor with A15 structure widely used in applica- 
tions. 

In the present paper we report experiments made on 
niobium nitride NbN films, another binary alloy shown 
here to have dendritic flux penetration in the supercon- 
ducting state. By combining MO imaging and magne- 
tometry we find threshold values for temperature and 
applied field, and a striking critical behavior in the size 
and morphology of the dendrites near the threshold. 

Thin films of NbN were fabricated by magnetron sput- 
tering on sapphire substrates. Two samples with thick- 
ness 0.16 and 0.29 pm were selected for the present stud- 
ies. Table 1 shows their critical temperatures T c , the 
width of the superconducting transition AT Cl the critical 
current density j c at 4.2 K, and the resistivity p^ in the 
normal state. All these parameters were obtained from 
transport measurements using a four-probe method. 

The magnetic moment was measured using a vibrat- 
ing sample magnetometer PARC with a He flow cryostat. 
Shown in Fig. 1 are magnetization curves m(H) for the 
NbN films obtained at 4.2 K. The instability manifests 
itself here in numerous and random jumps of m. The typ- 
ical jump amplitude is Am ~ 0.1m, which is much larger 
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TABLE I: Parameters of NbN films 
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FIG. 1: Magnetic moment as a function of applied field at 
4.2 K for two NbN films: 0.29 pm thick (upper) and 0.16 pm 
thick (lower). 



than the sensitivity of the magnetometer, 3d0 -5 emu. 
Magnification of a part of the m(H) curve (see the in- 
set) shows that the abrupt drops in magnetic moment 
are followed by a much slower increase before the next 
drop occurs. 

In the thinner film (lower panel) the jumps are seen to 
disappear above a threshold field H* as 1 kOe on the in- 
creasing field branch. Moreover, just below H* the jumps 
in m have a larger amplitude and occur less frequently 
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FIG. 2: Magneto-optical images of flux distribution in a 
0.29 nm thick NbN film at 3.5 K for increasing magnetic field. 
The image brightness represents the magnitude of the local 
flux density. 



than at low fields. On the descending field branch, the 
jumps reappear when m is fully reversed, and their am- 
plitude decreases rapidly. 

The nature of these jumps in the magnetic moment was 
clarified using MO imaging to visualize the dynamics of 
the full flux distribution^ The sample, with a Faraday 
active indicator film placed directly on top, was glued 
onto the cold finger of an optical cryostat, where it was 
cooled to 3.5-8 K in zero magnetic field (ZFC). Subse- 
quently, a perpendicular field was applied very slowly. 

For low fields most part of the superconductor is in the 
Meissner state and appear dark on the MO images. As 
the field increases, the flux penetrates gradually, start- 
ing preferentially from the weak places along the edges. 
At some field ss 10 Oe an abrupt invasion of a rela- 
tively large flux structure occurred, see Fig. [5J Further 
field increase resulted in formation of even larger and 
highly dendritic structures entering one by one. Eventu- 
ally, when reaching H w 28 Oe, the flux dendrites filled 
most of the film area. Upon further field increase, new 
dendrites continued to form, but now entering on top of 
already existing ones, as seen in the MO image taken at 
100 Oe. Note that the tilted lines of flux penetration 
near the edge is of a different origin, most likely due to 
polishing streaks in the substrate. 

Most features of the observed dendritic instability in 
the NbN films resemble those found previously in other 
materialsi 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i 1:L i 12 i 13 i 15 The dendrites propagate 
into the film faster than 1 ms, which is the time resolu- 
tion of the CCD camera recording our MO images. In 
fact, we expect that the propagation is even much faster, 
as was found using ultrafast MO imaging of dendrites in 
YBa2Cu307_,5 and MgB2 filmsAi^ Another characteris- 
tic feature is that once a dendritic structure is formed, it 
remains "frozen" and does not grow any further during 
subsequent increase of the applied field. Moreover, when 
the experiment is repeated under identical conditions, the 
exact field when the dendrites form and their exact pat- 
tern are never repeated. In the present films the dendritic 
instability was observed only below T* = 5.5 K, whereas 
a similar threshold temperature in MgB2 is 10 K£ Above 
these threshold temperatures the flux penetration is al- 




FIG. 3: First dendrites formed in a zero-field-cooled NbN film 
during increasing applied field H. Dendrites formed at higher 
temperatures are characterized by a stronger branching. 
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FIG. 4: Determination of fractal dimension of the first flux 
dendrite for different temperatures: 3.5, 3.8, 4.4 and 4.8 K 
from top to bottom. The curves are shifted vertically to avoid 
overlapping. 

ways spatially smooth and gradual in time. 

The instability disappears not only when T > T*, but 
also when the field becomes sufficiently high, H > H* , 
see Fig. flower panel). Our results clearly show that 
the threshold value depends on the field sweep direction. 
We propose that this dependence originates from vortex 
annihilation, which takes place only for the decreasing 
field case. Indeed, when H is increasing, the screening 
currents generate near the film edge a strong demagne- 
tization field of the same sign as H . However, for de- 
creasing H, the direction of screening currents and de- 
magnetization field changes to the opposite. As a result, 
the film edge experiences a negative external field that 
also partly penetrates inside. In this state, there is a line 
near the edge where vortices and antivortices meetiiii^ 
Their annihilation releases additional energy that can fa- 
cilitate the triggering of the instability^ Consequently, 
one may expect that the dendritic instability occurs in a 
wider range of applied fields H along the descending field 
branch as compared to the ascending branch. 

The existence of a threshold field H* was reported ear- 
lier in the magnetization studies of MgB2 filmsiSi Our ob- 
served asymmetry for the increasing and decreasing field 
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FIG. 5: Fractal dimension and total area of the first flux 
dendrite as functions of temperature. 

sweeps is also in agreement with results of Ref. 0, where 
dendritic jumps were found only for decreasing H . In- 
terestingly, during the dendrite growth the annihilation 
zone may propagate very deep into the film. This is con- 
firmed by observation of "negative" flux in the dendrite 
core that propagated into a film containing positive flux 
for a decreasing 

Figure [3] shows MO images of the very first dendrite 
formed in the 0.29 /im thick ZFC film during 4 experi- 
ments at slightly different temperatures. These dendrites 
were formed also at different (first jump) fields Hf 3 , and 
there is a clear tendency that flg increases with tempera- 
ture. Dendrites formed at higher T and Hf } - are also larger 
in size and more branching, a tendency can be traced all 
the way up to T* = 5.5 K. Note also that the m(H) curve 
in Fig. flower panel) exhibits increasingly larger jumps 
as the threshold field H* is approached. Therefore, it is a 
general trend that the dendritic structures have maximal 
size when the system is close to the stability limit, i.e. 
for H w H* or T w T*. 



To quantify these changes in morphology of the 
branching flux structures we made a fractal analysis of 
their shape. The MO images were discretized to obtain a 
cluster of pixels constituting a dendrite structure. Then 
we calculated the number of pixels N(R) that fall inside a 
circle of radius R with the center at the dendrite root. If 
the structure gives a power-law, N oc R D , the exponent 
gives the fractal dimension D of the cluster iSi 

The results of this analysis is presented in Fig.^J where 
the actual dendrite area represents N. We find a reason- 
ably good power-law behavior, and the fractal dimension 
increases with temperature, see Fig. [S] for a summary. 
The dimension changes from approximately unity at the 
lowest T to D = 1.75 for the most branching structure at 
4.8 K. Note that in general, the dendritic fingers do not 
exceed the middle of the strip, where the Meissner cur- 
rents (and hence Lorentz force) change to the opposite 
direction. This limits the size of all dendritic structures. 
As a result, more branching dendrites with larger D al- 
ways have a larger area. 

A similar temperature dependence of dendrite mor- 
phology has been reported earlier for Nb^i and for 
MgB^ samples. Various degrees of branching have also 
been obtained by simulations taking into account the 
heat produced by flux motion;2i2S which suggests a ther- 
mal origin of the instability. Our present results (i) give a 
quantitative measure of the branching, D, and (ii) show 
that simultaneous increase of D and dendrite area is a 
more general effect since it takes place when approach- 
ing the instability threshold H* (T) either by changing T 
or H. 

This work is supported by FUNMATOUiO, and The 
Research Council of Norway. 



E-mail for correspondence: t.h.johansen@fys.uio.no 

1 R. G. Mints and A. L. Rakhmanov, Rev. Mod. Phys. 53, 
551 (1981). 

2 A. VI. Gurevich, R. G. Mints and A. L. Rakhmanov, The 
Physics of Composite Superconductors, Begell House Inc., 
NY (1997). 

3 M. R. Wertheimer, J de G. Gilchrist, J. Phys. Chem Solids 
28, 2509 (1967). 

4 C. A. Duran, P. L. Gammel, R. E. Miller, D. J. Bishop, 
Phys. Rev. B 52, 75 (1995). 

5 M.S. Welling, R.J. Westerwaal, W. Lohstroh, R.J. Wijn- 
gaarden, Physica C 411, 11 (2004). 

6 P. Leiderer, J. Boneberg, P. Briill, V. Bujok, S. Herming- 
haus, Phys. Rev. Lett. 71, 2646 (1993). 

7 U. Bolz, B. Biehler, D. Schmidt, B.-U. Runge and P. Lei- 
derer, Europhysics Lett. 64, 517 (2003). 

8 T.H. Johansen et al, Europhys. Lett. 59, 599 (2002). 
|cond-mat/0104113| 

9 T.H. Johansen et al., Supercond. Sci. Technol. 14, 726 
(2001) . [condjmat/0108092] 

A.V. Bobyl et al AppL Phys. Lett. 80, 4588 (2002), 
cond-mat/0201260 



F. L. Barkov et a l., Phys. Rev. B 67, 064513 (2003). 
cond-mat/0205361 

S. C. Wimbush, B. Holzapfel, Ch. Jooss et al., J. Appl. 
Phys. 96, 3589 (2004). 

LA. Rudnev et al, Cryogenics 43, 663 (2003). 
|cond-mat/0211349| 

Ch. Jooss, J. Albrecht, H. Kuhn, S. Leonhardt and H. 

Kronmueller, Rep. Prog. Phys. 65, 651 (2002). 

For more publications on dendritic instability see 

http:/ /www. fys.uio.no/super/dend/ 

P. Leiderer, private communications. 

E. H. Brandt, and M. Indenbom, Phys. Rev. B 48, 12893 
(1993). 

E. Zeldov, J. R. Clem, M. McElfresh, and M. Darwin, 
Phys. Rev. B 49, 9802 (1994). 

M. R. Beasley, W. A. Fietz, R. W. Rollins, J. Silcox, and 

W. W. Webb, Phys. Rev. 137, A1205 (1965). 

Z. W. Zhao et al., Phys. Rev. B 65, 064512 (2002). 

J. Feder, Fractals (P lenum, New York, 1 988). 

I. S. Aranson et al. , |cond-mat /0407490 



